## Irr_q(m): number of monic irreducible polynomils
## of degree m over F_q

from PyM import *

def irr(q,m): 
    N = 0
    D = divisors(m)
    for d in D:
        N += mu_moebius(d)*q**(m//d)
    return N//m

  
M=[[irr(q,m) for m in range(1,8)] for q in [2,3,4,5,7,8,9]]


show(matrix(M))